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00 ■ ABSTRACT 



Context. We present S-band imaging of 18 low redshift (z < 0.3) BL Lac objects for which their host galaxies were previously 
resolved in the i?-band and the near-infrared //-band. For a subset of the objects, also U- and V-band imaging is presented. 
Aims. These multiwavelength data are used to investigate the blue-red-near-infrared colours and the colour gradients of the host 
galaxies of BL Lacs in comparison with other elliptical galaxies with and without nuclear activity. 
^ , Methods. For all the BL Lacs observed in the B- and (/-bands, and all objects at z < 0.15 in the f/-band, the host galaxy is clearly 

5_( . resolved. In all cases galaxies are well represented by an elliptical model, with average absolute magnitude Mb = -21.6 ± 0.7 and 

^ ' average scale length = 7.6±3.2 kpc. BL Lac host galaxies are therefore luminous (massive) elliptical galaxies in agreement of 

' previous studies in other bands. 

Results. The best-fit B-band Kormendy relation of (jig = 3.3 log (kpc) -I- 18.4 mag arcsec"^) is in reasonable agreement with 
that obtained for normal ellipticals and radio galaxies. This structural and dynamical similarity indicates that all massive elliptical 
galaxies can experience nuclear activity without significant perturbation of their global structure. The distributions of the integrated 
^ ' blue/near-infrared colour (with average B~H = 3.5+0.5) and colour gradient (with average A(B-R)/A(log r) = -0. 14±0.75) of the BL 

, Lac hosts are much wider than those for normal ellipticals, and most BL Lac objects have bluer hosts and/or steeper colour gradients 

00 . than those in normal ellipticals. 

' Conclusions. The blue colours are likely caused by a young stellar population component, and indicates a link between star formation 

04 ' caused by an interaction/merging event and the onset of the nuclear activity. This result is corroborated by stellar population modelling, 

indicating a presence of young/intermediate age populations in the majority of the sample, in agreement with low redshift quasar hosts. 
The lack of strong signs of interaction may require a significant time delay between the event with associated star formation episodes 
and the start of the nuclear activity. 

' Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: interactions - galaxies: nuclei - galaxies: photometry 

- galaxies: structure 
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H , 1. Introduction A nu mber of optical (e.g., 'Falomo & Kotilainen! 119991: 

■ - - lUrrv et al . 2000; Nilsson et al. 2003; Heidt et al. 2004) an d near- 

BL Lac objects are an extreme class of active galactic nuclei infrared (NIR) (e.g., Kotilainen et al. 1998; Scarpa etal. 20001 

(AGN) characterized by luminous, variable and polarized con- Cheung et al. 2003 ; Kotilai nen & Falomo 2004 ; Kotilai nen et al] 

tinuum emission across the electromagnetic spectrum and strong 2005|) imaging studies have shown that virtually all nearby (z < 

core-domi nated radio emission w ith apparent superluminal mo- 0.5) BL Lac objects are hosted in large and luminous elliptical 

tion (e.g.. iKollgaard et al.| [l992l). These properties suggest that galaxies, with average Mr 24.0 and average scale length 

they are strongly beamed objects dominated by synchrotron ^ 10 kpc. similarly to both FR I and FR II RGs dGovoni et all 
emission from a relativist ic jet aligned close to the line of sight |200Q) . BL Lac hosts are much brighter than L* galaxies (the 
dBIandford & Reeslll978l) . BL Lac objects have many similar- characteristic l uminosity of the Sch ec hter luminosity fi inctio n: 
ities to flat spectrum radio quasars (FSRQ) and they are often = -21.2, [Gardner et all (Il997h ; iNakamura etd] (|2003|) ). 
grouped together as blazars. According to the unified model While the morphology, such as jets and close companions, of 
of radio-loud AGN (UiTy & Padovani|[l995|), the parent popu- some of the hosts indicates a recent interaction, the large major- 
lation of BL Lac objects and FSRQs are low luminosity core- ity of them are indistii iguishable from in active massive ellipti- 
dominated FR I radio galaxies (RG), and high luminosity lobe- cals at similar redshift jScarpa et al.ll2000l) . 
dominated FR II RGs, respectively. Consequently, their orienta- 
tion independent properties, such as host galaxies and environ- Until recently, however, imaging of BL Lac host galaxies 
ments, should be identical to those of their parent populations. was obtained in one band only (usually T^-band). Therefore, 

little colour information exist for them, especially in the blue 

domain of the spectrum, where only a few objects have 

Send offprint requests to: T. Hyvonen been Studied in the Z?-band and even fewer in the f/-band. 
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iKotilainen & Falomol (|2004|) observed a sample of 23 low red- 
shift (z < 0.3) BL Lac objects in the //-band and combin ed with 
prev ious //-band imaging (Kotilainen et al. 1998; Scarp a et"aLl 
12000 ; Cheung et al. 2003) and optical /^-ban d data from litera- 
ture (.Falomo & Kotilainen, 19991 lUrrv et al]|200Q) . investigated 
the integrated R-H colours and colour gradients of a sample of 
41 BL Lac host galaxies. They found that BL Lac host galax- 
ies appear to be systematically bluer than inactive ellipticals and 
have a much wider distribution of host galaxy R-H colour (av- 
erage R-H - 2.2 + 0.4) and steeper colour gradient (average 
A(R-H)/A{logr) = -0.38 ± 0.60) than those for ina ctive elUpti- 
cal ga laxies with dominant old stellar population ( Peleti er et al.l 
Il990l) . Similar behaviour has been fo und for the colours and 
colour gradients in low redshift RGs dGovoni et al. I I2OOOI) and 
in low and intermediate redshift AGN and quasars with elhpti- 
cal ho sts (Schade et al. 2000; Orndahl et al. 2003; Jahn ke et al.] 
I2OO4I; ISanchez et al.l |2004). The blue colours are most likely 
caused by a young stellar population, indicating a recent star for- 
mation (SF) episode, possibly triggered by interaction or merg- 
ing between galaxies. The wide distribution of colours prob- 
ably reflects an object-to-object difference in the age of the 
most recent SF episode. Note that there is also spectroscopic 
evide nce for young/intermediate age popula tions in AGN hosts 
(e.g.. lNolan etaLll200ll;lRaimann et al.ll2005[) . 

In this study, we present multicolour optical imaging of a 
subsample of 18 BL Lac objects from the large, homogeneous 
sample of 41 sourc es for which high resolution R- and H- 
band imaging exists (jK otilainen et al.)ll998^, 'Scarp a et al.ll2000t 
ICheung et al.l2003l; |Kotilainen & Falomo 2004). Most of the ob- 
served objects have bluer R-H host colour than inactive ellipti- 
cals and blue band observations are paramount to assess whether 
their blue colours are caused by a young stellar population. All 
the 18 BL Lac objects were observed in the B-band, while for 
a subsample, we obtained also U- and V-band imaging. The U- 
band observations were mainly restricted to the most nearby tar- 
gets in the sample. Combining the deep high spatial resolution 
U-, B- and V-band imaging with the existing R- and //-band 
data, we are able to derive the colours and colour gradients over 
an extended wavelength range, where the //-band is sensitive to 
old stellar population that dominates the host galaxy mass, while 
the blue part of the spectrum traces the contribution from young 
stellar populations to the excess blue light in the host galaxies. 
We also derive structural properties, such as morphology and 
effective radius, as well as the absolute magnitude of the host 
galaxies in each observed band. The obtained UBVRH broad 
band colours are used in conjunction with stellar synthesis pop- 
ulation models to estimate the ages of the most recent SF episode 
in the host galaxies. 

In Section 2, we describe the sample, observations, data 
reduction, and methods of analysis. In Section 3, we present 
the results and discussion concerning the properties of the host 
galaxies. Summary and conclusions are given in Section 4. 
Throughout this paper, Hq = 70 kms 'Mpc^', Q,„ = 0.3 and 

- 0.7 cosmology is used. 



2. Observations, data reduction, and analysis 

The observations were carried out during several observing runs 
with different telescopes, with most of the observations done at 
tiie 2.5m Nordic Optical Telescope (NOT) with the ALFOS C in- 
strument using Bessel U, B and V broad band filters. Nine targets 
were observed only in the /J-band, five in the f/BV-bands and 
four in the t//?-bands. One target (MRK 421) was additionally 
observed in the //-band, to complement the sample presented in 



Table 1. Properties of the telescopes and instruments." 



Band Telescope+instrument Noi,j Scale 




arc sec px 


(1) (2) 


(3) (4) 


(/,S,V NOT/ALFOSC 


26 0.190 


U,B NOT/MOSCA 


3 0.217 


H NOT/NOTCam 


1 0.235 


U,B ESO NTT/EMMI 


3 0.273 



Column (1) gives the observed band; (2) the telescope and instru- 
ment used; (3) the number of observed targets, and (4) the spatial 
scale of the instrument. 



IKotilainen & Falomo (|2004. The B-band image of 3C 371, pub- 
lished in Nilsson et al. l (II997h . was kindly provided for us by K. 
Nilsson. A summary of the telescopes and instruments used is 
given in Table[T] Seeing during the observations varied between 
0.9 and 2.3 arcsec FWHM (average 1.5 + 0.5 arcsec and me- 
dian 1.3 arcsec). For each target, several integrations (typically 
of 600s duration) were obtained. The journal of the observations 
and properties of the sample are given in Table |2l 

Data reduction was performed in standard fashion using 
IRAFQ. This consisted of bias subtraction, flat field division and 
cosmic ray rejection. For each night, bias and flat field images 
were made from several median combined bias frames and twi- 
light flat field exposures and cosmic rays were rejected using 
IRAF procedures. Finally, individual images were aligned and 
combined to form the final image of the target. The observa- 
tions were mostly done in photometric conditions (see Table |2]i 
and several photometric standard stars from Landolt ( 1992) were 
observed during each night. Some targets were, however, ob- 
served during non-photometric nights and for these objects ad- 
ditional short exposures were subsequently obtained in photo- 
metric conditions to calibrate these frames using reference stars 
in the field. /T-correction from Poggianti (1997) was applied to 
the host galaxy magnitudes, but not to the nuclear magnitudes, 
since the nuclear component can be assumed to have a power- 
law spectrum (/v oc v^") with a ~ -\. Absolute magnitudes were 
also corrected for interstellar extinction calc ulated for e ach band 
from the /?-band extinction coefficient froml Urry et al.l (I2OOO). 

To derive the properties of the host galaxies, azimuthally av- 
eraged ID radial luminosity profiles were extracted for each BL 
Lac object and for a number of field stars. Any obvious extra fea- 
tures, such as nearby companions and/or foreground stars, were 
masked out from the image to avoid contamination of the radial 
luminosity profile. To obtain an accurate model for the nuclear 
region of the targets, it is important to have a well defined point 
spread function (PSF). In most cases there were a number of 
suitable stars in the relatively large field of view surrounding the 
objects to form a reliable PSF. The core and the wing of the PSF 
were derived from a faint and a bright star in the frame, respec- 
tively, and they were combined to form the final PSF model. This 
final PSF was compared with the profiles of individual stars in 
the frame to assure that this procedure resulted in a good and sta- 
ble representation of the true PSF. A representative case of this 
PSF comparison is shown in Fig.lTj where the individual stellar 
profiles are in good agreement with the adopted PSF model well 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Name 




V 


Mb 


Filter 

(sec) 


eip 


Date 
(arcsec) 


FWHM 


Photometric? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


lES 0229+200 


0.139 


18.0 


-21.7 


U 


2700 


12/11/2004 


1.7 


Y 










B 


900 


12/11/2004 


1.9 


Y 










V 


1800 


12/11/2004 


2.0 


Y 


PKS 0521-365 


0.055 


14.6 


-22.3 


B 


480 


14/03/2005 


1.4 


Y 


PKS 0548-322 


0.069 


15.5 


-22.0 


B 


900 


15/03/2005 


1.3 


Y 


IH 0706+591 


0.125 


19.5 


-21.0 


U 


1600 


27/01/2004 


1.2 


Y 










B 


1500 


27/01/2004 


0.9 


Y 










V 


600 


27/01/2004 


1.2 


Y 


MRK 421 


0.031 


13.8 


-22.9 


U 


900 


08/03/2005 


0.9 


N 










B 


300 


08/03/2005 


0.9 


N 










H 


150 


25/03/2005 


1.8 


Y 


MRK 180 


0.045 


15.0 


-22.1 


U 


900 


06/05/2002 


1.2 


Y 










B 


1600 


10/03/2003 


1.2 


Y 


lES 1212+078 


0.136 


16.1 


-23.6 


B 


1800 


08/03/2005 


0.9 


N 


MS 1229.2+6430 


0.164 


16.9 


-21.7 


B 


1800 


10/03/2004 


2.1 


Y 


lES 1255+244 


0.140 


15.4 


-24.3 


B 


1800 


09/03/2005 


1.9 


Y 


PG 1418+546 


0.152 


15.7 


-23.7 


U 


900 


07/05/2002 


2.1 


Y 










B 


700 


07/05/2002 


1.3 


Y 


lES 1426+428 


0.129 


16.5 


-22.5 


B 


1600 


10/03/2003 


1.8 


Y 


lES 1440+122 


0.162 


17.0 


-22.9 


B 


1800 


08/03/2005 


0.9 


N 


AP LIBRAE 


0.049 


14.8 


-21.7 


U 


700 


07/05/2002 


2.0 


Y 










B 


900 


15/03/2005 


1.1 


Y 


MS 1552.1+2020 


0.273 


17.7 


-23.5 


B 


1200 


09/03/2005 


1.9 


Y 


MRK 501 


0.034 


13.8 


-22.4 


U 


600 


07/05/2002 


1.2 


Y 










B 


240 


08/03/2005 


1.1 


N 










V 


loU 


Uo/Uj/ZUUJ 


1 9 
1 .z 


NT 


IZw 187 


0.055 


16.4 


-21.1 


U 


500 


07/05/2002 


2.0 


Y 










B 


900 


25/09/2003 


1.1 


Y 










V 


300 


08/03/2005 


1.3 


N 


3C 371 


0.051 


14.2 


-22.6 


U 


1200 


01/10/2000 


2.1 


Y 










B 


3600 


02/05/1997 


0.8 


Y 










V 


480 


08/03/2005 


1.3 


N 


BL LACERTAE 


0.069 


14.7 


-22.4 


B 


720 


01/10/2000 


2.3 


Y 



" Column (1) gives the name of the BL Lac object; (2) the redshift; (3) the V-band apparent magnitude; (4) the S-band absolute magnitude; (5) 
the filter used; (6) the total integration time; (7) the date of the observation; (8) the seeing FWHM and (9) photometric (Y) or non-photometric 
(N) conditions. 



into the domain where the host galaxy becomes dominant over 
the PSF. The only exception was the field of MRK 42 1 where 
there are no suitable stars available. In this case, the PSF was 
estimated using standard stars observed during the same night 
with similar seeing conditions. 

The luminosity profiles were decomposed into a point source 
(represented by the PSF) and an elliptical galaxy components by 
an iterative least-squares fit to the observed profile. There are 
three free parameters in the fit: the PSF normalization, the host 
galaxy normalization and the effective radius of the host galaxy. 
The data were fit using the r''''* de Vaucouleurs law for ellipti- 
cal galaxies to represent the host galaxy. The host galaxy was 
considered to be resolved if the PSF + host galaxy fit resulted in 
a considerably lower value than the PSF fit only. The uncer- 
tainty in the derived host galaxy magnitudes was estimated to be 
~ 0.2 mag. 



3. Results 

From the images ID azimuthally averaged radial luminosity pro- 
files of the BL Lac objects were extracted in the U-, B- and V- 



band. The profiles together with the best-fit model overlaid are 
presented in Appendix in online version of the paper 

We were able to clearly resolve the host galaxy in all objects 
in the B- and V-bands, and in 8/9 objects in the [/-band. The only 
target that remained unresolved in the [/-band, PG 1418+546, is 
one of the most distant objects in the sample. The results of the 
modelling are summarized in Table [3] 

3.1. Luminosities and sizes of fhe hosf galaxies 

Fig.|2]shows the distribution of the B-band absolute magnitudes 
of the BL Lac host galaxies (this work), low redshi ft radio-loud 



RLQ) and radio-quiet (RQQ) q uasar hosts from Jahnke et al.l 
I2OO4I) . low redshift RGs from Govoni et al.' (200(f) an d inac- 
tive el liptic al galaxies fr om Peletieret al. (1990), Colber tet al.l 
(I2OOII) and iBower et all ( 119921) . The average B-band absolute 
magnitudes of the BL Lac objects, quasar and RG host galax- 
ies in these samples are presented in Table H] where all absolute 
quantities were transformed into the cosmology adopted here. 

Since the selection criteria of the various samples are some- 
what non-homogeneous, there is a possibility of selection ef- 
fects. They are likely to affect absolute magnitudes but since 
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Table 3. Properties of the host galaxies.^ 



Name 


Filter 


A 


z 












jut 


A/I 


hi 1 u 






mag 




mag 


mag 


mag arcsec~ 


arcsec 


kpc 


mag 


mag 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


lES 0229+200 


U 


0.84 


0.139 


19.22 


18.83 


20.94 


2.75 


6.55 


-20.6 


-21.5 


0.44 




B 


0.71 




19.58 


18.57 


22.39 


5.65 


13.4 


-20.2 


-21.9 


0.21 




V 


0.54 




19.67 


17.58 


21.72 


4.90 


11.7 


-19.9 


-22.2 


0.12 


PKS 0521-365 


B 


0.25 


0.055 


16.58 


16.17 


21.50 


7.60 


8.02 


-20.6 


-21.3 


0.52 


PKS 0548-322 


B 


0.17 


0.069 


17.42 


16.11 


22.31 


9.40 


12.2 


-20.2 


-21.8 


0.23 


IH 0706+591 


U 


0.50 


0.125 


17.67 


18.55 


21.90 


3.60 


7.83 


-21.6 


-21.1 


1.6 




B 


0.42 




18.38 


17.86 


21.21 


3.50 


7.63 


-20.8 


-22.0 


0.33 




V 


0.32 




18.10 


16.57 


20.40 


3.65 


7.96 


-21.0 


-22.8 


0.19 


MRK 421 


u 


0.16 


0.031 


13.41 


15.22 


21.29 


7.80 


4.80 


-22.4 


-20.7 


4.8 




B 


0.14 




14.92 


15.43 


21.78 


8.90 


5.48 


-20.9 


-20.5 


1.4 




H 


0.00 




15.18 


14.29 


16.82 


6.95 


4.45 


-20.5 


-24.9 


0.02 


MRK 180 


U 


0.16 


0.045 


14.83 


15.81 


20.83 


4.95 


4.33 


-21.8 


-21.0 


2.1 




B 


0.14 




16.44 


15.45 


21.91 


9.45 


8.27 


-20.2 


-21.4 


0.33 


lES 1212+078 


B 


0.13 


0.136 


19.17 


18.53 


22.09 


3.80 


8.89 


-19.9 


-21.3 


0.28 


MS 1229.2+6430 


B 


0.17 


0.164 


19.14 


18.19 


21.28 


3.15 


8.57 


-20.4 


-22.2 


0.19 


lES 1255+244 


B 


0.13 


0.140 


19.45 


18.33 


21.90 


3.55 


8.50 


-19.7 


-21.6 


0.17 


PG 1418+546 


U 


0.09 


0.152 


15.95 










-23.4 








B 


0.08 




17.20 


18.14 


20.25 


1.90 


4.87 


-22.1 


-21.9 


1.2 


lES 1426+428 


B 


0.13 


0.129 


17.67 


17.50 


19.35 


1.40 


3.13 


-21.3 


-22.1 


0.48 


lES 1440+122 


B 


0.13 


0.162 


18.14 


19.17 


22.43 


3.65 


9.84 


-21.4 


-21.2 


1.2 


AP LIBRAE 


U 


0.84 


0.049 


15.26 


16.10 


19.07 


2.55 


2.42 


-22.2 


-21.6 


1.7 




B 


0.71 




16.88 


16.52 


19.74 


2.90 


2.75 


-20.5 


-21.1 


0.58 


MS 1552.1+2020 


B 


0.29 


0.273 


19.22 


19.07 


20.09 


3.40 


13.53 


-21.7 


-23.1 


0.28 


MRK 501 


U 


0.16 


0.034 


14.59 


14.43 


21.07 


8.90 


5.98 


-21.4 


-21.7 


0.76 




B 


0.13 




14.89 


14.20 


20.76 


9.45 


6.35 


-21.1 


-22.0 


0.44 




V 


0.10 




14.43 


13.00 


19.68 


9.45 


6.35 


-21.5 


-23.2 


0.21 


IZw 187 


U 


0.25 


0.055 


16.53 


17.56 


21.11 


2.65 


2.80 


-20.5 


-19.8 


1.9 




B 


0.21 




17.12 


17.14 


21.65 


4.35 


4.59 


-20.0 


-20.3 


0.76 




V 


0.16 




16.39 


15.82 


20.93 


5.15 


5.43 


-20.7 


-21.4 


0.52 


3C 371 


u 


0.45 


0.051 


15.33 


16.17 


22.13 


8.95 


8.80 


-21.9 


-21.2 


1.9 




B 


0.38 




15.86 


15.98 


21.12 


6.15 


6.05 


-21.3 


-21.4 


0.91 




V 


0.29 




14.98 


14.94 


20.12 


5.60 


5.51 


-22.1 


-22.2 


0.91 


BL LACERTAE 


B 


1.61 


0.069 


14.71 


17.35 


19.79 


3.20 


4.15 


-24.3 


-22.0 


8.32 



" Column (1) gives the name of the object; (2) the filter; (3) the interstellar extinction in the U, B- and V-bands; (4) the redshift of the object; (5) 
the apparent nuclear magnitude; (6) the apparent host galaxy magnitude; (7) surface brightness at the effective radius (8) and (9) the apparent 
and absolute effective radius, respectively; (10) the absolute nuclear magnitude; (11) the absolute host galaxy magnitude and (12) the absolute 
nuclear/host luminosity ratio. 



the main focus of this study is on the colour properties of BL 
Lac hosts as compared with other ellipticals, our main results 
are only marginally affected by selection effects. Furthermore, 
our BL Lac sample is well matched to RGs and quasars in red- 
shift (see Table |4]l to minimize any bias introduced by possible 
evolutionary effects. 

The average B-band absolute magnitude of the 18 BL Lac 
hosts is Mb - -21.6 + 0.7, i.e. ~1 mag brighter than for low 
redsh ift quasar hosts (average Mb - -20.9 + 0.4 jjahnke et al.l 
12004 ). All the BL Lac hosts are in the luminosity range between 
M; + 0.5 andM!, + 2.5 where M* = -20.8 (Gardner et al. 199^ 
iNakamura et al.ll2003l) is the characteristic luminosity of the 
Schechter luminosity function for elliptical galaxies. The aver- 
age luminosity of the BL Lac hosts is in agreement with inactive 
isolat ed early-type galaxies (Mb - -21.5 + 0.7; (Colb ert et al.l 

200 Ih) but is ~ 1 '^^S brighter than ellipticals observed by 

Peletier et all d 19901) (Mb = -20.9 + 1.6). The distribution of 
inactive ellipticals has a tail toward fainter magnitudes that is 
not present in the BL Lac sample. 

According to the unified model of AGN dUrrv & Padovanil 

119951) . BL Lac objects are physically similar to FR I RGs but 



viewed from a different orientation. The classification of BL Lac 
objects with FR Is is based on their radio morphology, show- 
ing that BL La c objects have similar exte nded radio properties 
to FR Is (e.g.. lAntonucci & Ulvestadlll985 ) but, however, there 
are some indications that the radio properties of at least some 
BL Lac objects are similar to FR lis rather than FR Is (e.g., 
Kollgaar d et"aLlll992l:IStanghellini et anil997h . There is also ev- 
idence that some BL Lac object s share the luminosity properties 
of FR lis (ICassaro et al. l [T999h . The average B-band luminos- 
ity of low redshift FR I and FR II RGs is Mfri ^ -22.1 ± 0.6 
and M pRii - -2 1 .9 + 0.6. respectively, obtained bv lGovoni et al.l 
(2000). BL Lac hosts appear, therefore, to be slightly fainter than 
both FR Is and FR lis, but based on Kolmogorov-Smirnov statis- 
tics, the luminosity distributions for the BL Lac hosts and the 
combined sample of RGs are indistinguishable. Our result sug- 
gests that both FR I and FR II RGs can be considered as a parent 
populati on of BL L ac objects, consistent with the result obtained 
bv palomo & Kotilainen (1999I) . 

The average effective radius of the BL Lac host galaxies in 
the B-band is = 7 .6 ± 3.2 kpc, consistent with results from 
previous studies (e.g.. lNilsson et al]|2003l : iFalomo & KotilainenI 
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Table 4. Comparison of the average fi-band host galaxy properties.^ 
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Name A' < z> < Mgihost) > < > 

mag kpc 

(1) (2) (3) (4) (5) 



BL Lacs (this work) 


18 


0.104 


+ 


0.064 


-21.6: 


tO.7 


7.6 4 


3.2 


Quasars (Jahnke et al. 2004) 


12 


0.094 


+ 


0.043 


-20.9 : 


tO.7 


6.0 4 


-2.7 


FR I RGs (Govoni et al. 2000) 


17 


0.069 


+ 


0.031 


-22.1 : 


1:0.6 


15.3 


t9.2 


FR II RGs (Govoni et al. 2000) 


6 


0.050 


+ 


0.016 


-21.9: 


1:0.6 


11.4 


t7.7 



^ Column (1) gives the sample; (2) the number of objects; (3) the average redshift; (4) the average host galaxy magnitude, and (5) the average 
effective radius (< > for RGs are in i?-band). 




0.5 ^ ;,:\, :■'■[ r,^ ^ 

-0.5 ^ i"'"'';:' 'i ^ 

5 10 15 

Radius (arcsec) 

Fig. 1. Upper panel: Comparison of the PSF profile extracted 
from field stars (the solid line) with the profiles of individual 
stars (dotted lines) in the field of BL Lac object. Lower panel: 
The difference between the profiles of individual stars and the 
PSF model. 



1 19991) that BL Lac objects are hosted in large elliptical galax- 
ies. BL Lac hosts are of similar size in the fi -band to quasar 
host galaxies (Re - 6.0+2.7 kpc; iJahnke et all ^2004,) ) but they 
are shghtly smaller than RGs in /J-band (Re = 13.6 ± 8.7 kpc); 
iGovoni et al.l (l2000h ). Note, however, that the latter comparison 
may be biased by a luminosity-dependent selection effect. 

It is well known that the effective radius of inactive 
late-type galaxies increases toward shorter wavelengths (e.g.. 



iMollenhoff & Heid3l200ll) . This dependence delineates the mor- 
phology of SF regions in the galaxies, in the sense that SF (blue 
bands) occurs in the spiral arms while the old stellar population 
(red bands) dominates the bulge component. On the other hand, 
|de Griis ( 1998) found no change in the effective radius from B- 
to /T-band for inactive early-type galaxies (up to Sa spirals). The 
average effective radii of our BL Lac host galaxy sample are 
Re(U) = 6. 4+2 3 kpc, Re(B) = 7.6+3.4 kpc and/?,(y) = 7.4+2.6 
kpc, comparable to those in the //-band (Re(H) = 7.8 + 4.3 kpc; 
iKotilainen & Falomo 2004)) and in the /?-band (Re(R) = 9 + 5 
kpc; (iFalomo & Kotilainen 1999)). In agreement with the result 
for inactive ellipticals de Griis, (,1998i) . we find that the effective 
radius of the early-type AGN host galaxies does not decrease 
with wavelength. If anything, one can note a slight increasing 
trend. However, note that this result should be taken with cau- 
tion because of the relatively large errors in determining the ef- 
fective radius depending on the well known degeneracy between 
surface brightness and effective radius. 

3.2. Kormendy relation 

According to the Kormendy relation, there is a tight relation 
between the effec t ive radius Re and the surface b rightness /i^ 
(iKormendvl [T977t iKormendv & Diorgovskil IT989h . This rela- 
tion is a 2D projectio n of the 3D F undamental Plane (e.g., 
iDressler etafll9 87; Diorgovski & Davi s 1987) that Unks Re and 
fie with the stellar velocity dispersion cr. This relation, well es- 
tablished for inactive elliptical ga laxies and the bulges of spira l 
galaxies in nearby clusters (e.g.. iJorgensen et"aLlll993L Il996l) . 
is related to the morphology and the dynamical structure of the 
galaxies and gives important information about their formation 
processes, indicating underlying regularity within the galaxy 
populations. 

Fig. [3] presents the B-band Kormendy relation for our sam- 
ple o f 18 BL Lac obje cts, compared with that for nearby 
RG s dGovoni et all l2000 l). isolated inactive elliptical galaxies 
(Re da et al. Il2005h and inactive ellipticals in the Coma cluster 
(Jorgens en et alJll993h . The corresponding surface brightnesses 
jUp and effective radii Rg for the BL Lac hosts are given in Table[3] 
The surface brightnesses of the host galaxies were corrected for 
Galactic extinction and cosmological dimming (lOx log(l+z)). 
Note that the surface brightnesses fie and effective radii Re refer 
to the isophote that contains half of the total luminosity of the 
galaxy. Another commonly used definition for fie and r^. is to de- 
rive these values directly from the de Vaucouleurs fitting of the 
galaxy, but that definition is model dependent. 

The best-fit linear relation for the BL Lac hosts is fie = 3.3 
log Re (kpc) + 18.4 mag arcsec"^. It is somewhat steeper than 
that for inactive elliptical galaxies and RGs. It is, however, con- 
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Fig. 3. The B-band He-Re Kormendy relation for the BL Lac 
host galax ies (filled symbo l s), and for i nactive elliptical g alaxies 
(asterisks; iJorgensen et"!!] 11993') and lReda etall (|2005)). The 
solid, dotted and long dashed lines represent the best line ar fit 
for BL Lac hosts, inactive ellipticals jJorgensen et alj ( 19931) an d 
iReda et alj (l2005h ). and low redshift RGs jFasano et alJ (Il996l) ). 
respectively. 



rienced a phase of nuclear activity in the past with little influence 
on the global structure of the galaxy. 



Fig. 2. Histogram of the absolute B-band magnitudes of the BL 
Lac hosts (top panel; this work), low redshift {z < 0.2) RLQ 
and RQQ hosts (second panel; (Jahnke et al. 2004)), low red- 
shift RGs (third panel; (Govoni et al. 2OO( )0) and inactive ellip- 
tical g alaxi es (bottom panel; | Peletier et al.lll990l) ; IColbert et al.l 
(1200 11) and lBower et all (Il992h ). The vertical short dashed line 
represent the luminosity of the L* galaxies. 



sist ent with the relation ob t ained bv lFalomo & KotilainenI d 19991) 
and lKotilainen & Falomol (|2004|) for BL Lac hosts in the R- and 
i/-bands, indicating that the dynamical structure of the hosts 
does not change with wavelength. Note that one BL Lac object 
(MS 1552.1H-2020, filled triangle in Fig. IS clearly deviates from 
the He - .Rf-relation of the sample, towards large effective radius 
and high surface brightness. Note that MS 1552. 1 -1-2020 also has 
one of the bluest host galaxies in the sample (B-H-3.2; Table|5]i, 
and its displacement from the Kormendy relation may indicate 
recent S F as suggested fo r similar cases among inactive ellip- 
ticals bv lReda et al.l (|2005|) . Overall, it seems that BL Lac hosts 
are not very far from the Fundamental Plane of inactive elliptical 
galaxies, indicating that they are dynamically similar It supports 
the view that BL Lac hosts and normal ellipticals have a similar 
formation process and that all elliptical galaxies may have expe- 



3.3. Host galaxy colours 

The host galaxies of all the BL Lac objects in our sample were 
previously studied in the opt i cal j^-band and in the N IR //-band 
(iFalomo & KotilainenI [19991 iKotilainen & Falomol [2004 ). The 
new observations presented here represents the first host galaxy 
study performed in the blue part of the spectrum. The combina- 
tion of blue, optical and NIR data allows us to assess the issue of 
the optical-NIR colours of the BL Lac host galaxies. While the 
//-band is sensitive to old stellar populations, the blue domain 
of the spectrum is especially important in studying SF and the 
presence of a young stellar population component in the galax- 
ies. 

The integrated rest-frame U-B, B-V and B-R colours of 
our BL Lac sample are given in Table [5| The average colours 
of the BL Lac hosts, compared to those obse rved in RLQ 
and RQQ host galaxies ( Jah iike et aT] |2004). RGs ( Govoni et al] 



200 0), and inactive elliptical galaxie s CPeletier et al. 119901; 



^ il ga 

Bo wer et all [19921; [Colbert et all [200 lb. and flieoretical l y pre- 
dict ed for inactive e llipticals ([Fukugita & Ichikawd [l995[; 
Fioc & Rocca-Volmera nge 1999) are presented in Table [6| 
Because [/-band data is available only for eight objects in our 
sample, we prefer to use the B-H colour as the longest baseline 
colour 

The average B-H colour of the BL Lac hosts is B-H=3.5 + 
0.5 that is sUg hdy redd er than that of RLQ and RQQ hosts (B- 
H^2.9 + 0.3, Jahnke et a l. (2004)) but slig htly bluer flia i i that of 
inactive elliptical galaxies (B-H ^3 .8+0.3 . iColbert et al.[ ([200 1[) ) . 
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Fig. 4. The distribution of the B-R colour for BL Lac hosts (top 
panel; this work), low redshi ft (z < 0.2) RLQ and RQQ hosts 
(seco nd panel; [ja hnke et alj (|2004 ). low redshift RGs (third 
panel; iGovoni et al. (2000)) and inact ive elliptical galaxi es (bot- 
tom panel; , Colbert et al.. (.2001.) and Peletier et alj dl990h ). 



The same situation applies when considering the B-V and B- 
R colours (Table |6]l. In fact, the optical colours of the BL Lac 
hosts are very similar to those of intermediate/late-type (Sb - 
Sbc) inactive galaxies with significant ongoing SF (Table |6l). 
Fig.|4]shows the distribution of the B-R colour for BL Lac hosts, 
quasar hosts (Jahnke et al. 2004.). RGs (iGovoni et al., 2000) and 
inactive ellipticals ( IColbert et al.l200ll) . BL Lac hosts clearly ex- 
hibit a much wider colour distribution and have a significantly 
bluer host galaxy population that is not present for inactive el- 
lipticals. 

It has long been known that the integrated colours of ellip- 
tical galaxies become redder toward higher lum inosity (mass). 
(e.g.. lBower et all 1 9921: iKodama & Arimotoll997ft . This colour- 
magnitude relation links the properties of the stellar populations 
of early-type galaxies with their structural properties and pro- 
vides important information about their formation and evolu- 



tion. It depends on the combined effects of age and metallicity 
on the dominant stellar population (more massive galaxies are 
both older and more metal-rich than less massive galaxies). It is 
identical in different galaxy clusters, such as in Virgo and Coma 
clusters dBower et a l.^^1992") and overall shows little dependence 



on environment TTerievich^ etal] 120011; [B ernardi et al. 2003] 



and redshift (e.g., Aragon-Salamanca et alJll993i:iKodama et akl 
Q998; Holden et al. 2004), indicating that massive ellipticals 
formed in an intense starburst at high redshift followed by pas- 
sive evolution, with no major SF episodes since z ~2. 

The B-R vs. R and B-H vs. H colour-magnitude diagrams 
for th e BL La c hosts (this work), low redshift RLQ and RQQ 
hosts dJahnke et al. 2004), R Gs (Govoni et al. 2000) , inactive el- 
lipticals (Peletier et al. 1990; IColbert et alJl200lT and ellipticals 
in Virgo and Coma clusters dBower et al.ll 1992 b are presented 
in Fig. |5] It is evident that the BL Lac hosts do not follow 
the relatively tight colour-magnitude relation of inactive ellip- 
ticals. Instead, they have a significantly broader colour distribu- 
tion and the majority of them appear to be bluer than inactive 
ellipticals of similar luminosity. Indeed, such colours are more 
similar to those found in intermediate/late-type inactive galax- 
ies that have significant recent S F. This result is consistent with 
the co lours of low redshif t RGs (IGovoni et alJl200Q) and quasar 
hosts dJahnke et al.ll2004l) . Although this colour difference may 
indicate that the colour-magnitude relation for elliptical galaxies 
breaks do wn at high luminosities as sug gested by Govoni et al] 
d2000l) and lKotilainen & Falomol d2004l) . note that the elliptical 
galaxies cover the same luminosity range as the BL Lac hosts. 
A more likely explanation for the blue colours of BL Lac hosts 
(and the other AGN hosts in the diagrams) is that they have expe- 
rienced recent SF. The wide colour distribution indicates a range 
of timescale since the latest SF episode, such that the bluest hosts 
have experienced the most recent SF whereas the reddest hosts 
have experienced little or no SF and are dominated by an old stel- 
lar population similar to those in inactive ellipticals. Especially 
note that there is one BL Lac host exhibiting red B-R colour (B- 
R > 2.0) and two hosts having red B-H colour (B-H > 4.4), red- 
der than any of the inactive ellipticals. From this point of view, 
the bluest hosts may hold evidence from an event that triggered 
both the nuclear activity and the strong SF. This is consistent 
with the blue coloii rs (B-R = 1.1 + 0.1) in RLQ and RQQ hosts 
(iJahnke et al.ll2004l) . Simi larly blue host colou rs were also found 
in low redshift quasars bv lSchade et"aLld2000t) . 



3.4. Colour gradients of the host galaxies 

It is well known that nearby inactive elliptical galaxies do 
not have uniform colour but instead have negative colour gra- 
dients, i.e. they b ecome bluer with increasing radius (e.g., 
IPeletier et al.lll990h . These colour gradients have been widely 
interpreted as due to radial variations in the du st content and/or 
SF history of the galaxies (e.g., [Kim 1989; Pel etier et al]| 19901; 
iGoudfrooii et alJll9 94; Michard 2000). As dust reddening is rel- 
atively s mall in elhpticals , while normal early-type galaxies can 
be dusty dTran et alj|200ll) . most of the colour gradients are usu- 
ally ascribed to the combined effect of stellar metallicity and 
age gradients across the galaxies, with the outer regions being 
younger and/or having lower metallicity. 

For the BL Lac host galaxies, we derived the radial U- 
B, B-R and B-H colour gradients u sing the R- and //-ban d 
host galaxy luminosity profiles from Kotilainen et al. (1998!); 
Fal omo & K o tilaine n (1999); Kotilainen & Falomo (2004^; 
iScarpa et al.l (I2OOOI) . The B-R colour profiles are presented in 
Fig. |6] and all colour gradients are presented in Table |5] Each 
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Table 5. Colours and colour gradients of the BL Lac host galaxies." 



Name U-B B-V B-R R~H A(U-B)/A(logr) K{B-R) I MJogr) A(B~H)/ A(logr) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


lES 0229+200 


0.4 


0.3 


2.0 


2.5 


0.92 


-0.44 


-0.17 


PKS 0521-365 




0.9 


0.9 


2.3 




-1.02 


-1.66 


PKS 0548-322 






1.1 


2.1 




1.10 


-0.90 


IH 0706+591 


0.9 


0.8 


1.3 


2.2 


-0.01 


0.17 


-0.12 


MRK 421 


-0.2 




2.2 


2.2 


0.17 


-0.91 


-0.31 


MRK 180 


0.4 




0.5 


2.7 


0.75 


-1.24 


-0.91 


lES 1212+078 






1.9 


2.2 




0.48 


-0.85 


MS 1229.2+6430 






1.3 


0.9 




-0.12 


-0.74 


lES 1255+244 






0.9 


2.1 




-1.08 


-0.91 


PG 1418+546 






1.4 


2.2 




-0.21 


0.86 


lES 1426+428 






0.8 


2.5 




0.90 


1.19 


lES 1440+122 






1.6 


2.4 




0.16 


-0.11 


AP LIBRAE 


-0.5 




1.1 


2.2 


0.17 


0.43 


-0.31 


MS 1552.1+2020 






0.5 


2.7 




-0.41 


-0.37 


MRK 501 


0.4 


1.2 


1.1 


2.4 


0.01 


0.64 


-0.62 


IZw 187 


0.5 


1.1 


1.2 


2.3 


0.63 


-0.32 


-0.54 


3C 371 


0.2 


0.8 


1.5 


1.9 


-0.44 


-1.26 


-1.51 


BL LACERTAE 






1.0 


2.4 




0.63 


-1.70 



" Column (1) gives the nam e of the object; (2) the U-B colour; (3) the B-V colour; (4) the B-R colour; (5) the R-H colour (from 
iKotilainen & Falomol i2004h ^; (6) the U-B colour gradient; (7) the B-R colour gradient and (8) the B-H colour gradient of the host galaxy. 

Table 6. Average colours of the host galaxies.'' 

The sample U-B B-V B-R R-H B-H 

(1) (2) (3) (4) (5) (6) 



BL Lacs (this work) 0.3 ± 0.4 0.8 ± 0.2 1.2 ± 0.5 2.2 ± 0.4 3.5 ± 0.5 

Quasars f Jahnke et al. 2004) 0.5 + 0.1 1.1 +0.1 1.9 + 0.2 2.9 + 0.3 

RGs (Govoni et al. 2000) 1.4 ± 0.2 

Ellipticals (Peletier et al. 1990) 0.6 ± 0. 1 1.6 ± 0. 1 

Ellipticals (Bower et al. 1992) 3.8 ± 0.2 

Ellipticals (Colbert et al. 2001) 1.5 ± 0.2 

Ellipticals (Fukugita & Ichikaw a 1995) 1.0 1.6 

Ellipticals (Fioc & Rocca-Volmerange 199^ 2.2 3.8 

Spirals (Sb) (Fukugita & Ichikawa 1995; Fioc & Rocca-Volmerange 1999) 0.7 1.2 2.1 3.3 



^ Column (1) gives the sample; and columns (2) - (6) the average U-B, B-V, B-R, R-H and B-H colour, respectively. 



gradient was cut at the radius < 2 arcsec to avoid any contam- 
ination of the host galaxy from the bright nucleus. The colour 
as a function of logarithmic radius is usually well represented 
by a linear fit. As expected for early-type galaxies, most of 
the BL Lac host galaxies exhibit a negative colour gradient 
with average A{B-R)/ A{logr) = -0.14 ± 0.75. The same neg- 
ative trend was also found for the R-H c olour gradients of BL 
Lacs ( A(R-H)/A(lof^r) = -0.38 + 0.60 bv lKotilainen & Falomol 
( |2004 ). The amplitude of the B-R colour gradients of the BL 
Lac hosts is consistent with the average B-R gradient of RGs 
MB-R)IMlogr) = -0.16 + 0.17 jGovoni et al.i.2000) and inac- 
tive ellipticals A(B-R) / A(logr) = -0.09 + 0.07 jPeletier et alj 
11990 ). but the distribution is much wider than that of the RGs 
and inactive ellipticals. Similar distribution was also obtained 
for the R-H gradient of BL Lac hosts by iKotilainen & Falomol 
(12004 . 

However, note that some BL Lac host galaxies have little 
evidence for colour variation or even have an inverted (posi- 
tive) colour gradient, the steepest of them in PKS 0548-322, 



lES 1426+428 and MRK 501. Similarly steep positive R-H 
colour gradient for lES 1426 +428 was previously observed by 
IKotilainen & Falomol ( l2004i) . Since the B-H colour is more 
sensitive to the dust content than the R-H colour, it indicates 
that the positive B-R and R-H gradients of this target are 
due to radial variations in the dust content of its host galaxy. 
Indeed, recent HST observ ations of RGs (e.g., iMartel et aTl 
I2OOOI : IrVemblav et alj|2007h have shown that there is a signif- 
icant number of RGs that contain dust in a variety of spatial 
distributions, such as circumnuclear disks and dust lanes at kpc 
scales. On the other hand, dust is not a reasonable explanation 
for e.g. MRK 501 which has a positive B-R but a negative R- 
H gradient. In this case, the inverted profiles indicates SF in the 
inner region of the host galaxy because the host galaxy is well re- 
solved. That indication is also supported by recent spectroscopic 
observations of ongoin g SF in the nuclear r egions of a BL Lac 
object (PKS 2005-489: IBressan et al.1 (120061) ). 

Fig. [T] shows the U-B colour profiles for the eight BL Lac 
host galaxies for which we have the available data (Table |5]l. 
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Fig. 6. The B-R co lour profiles for the BL Lac hos t galaxies, derived using the 7? -band data from iFalomo & KotilainenI (Il999t) : 
IScarpa et alj (|2000|) and lKotilainen & Falomol(l2004 . 



All the host galaxies, except one (3C 371), have a positive 
colour gradient. 3C 371 has some signature of recent interaction 
(iNilsson et alJI 19971) which might be the reason for the negative 
U-B colour gradient. Average colour gradient for all eight ob- 
jects is A{U-B)/A{Iogr) - 0.28 + 0.46, although for four objects 
the gradient is very flat. 



3.5. Stellar population model fits 

The UBV colours combined with the previously observed R- and 
//-band data can be used in conjunction with stellar synthesis 
population models to estimate the approximate ages of the recent 
SF episodes in the BL Lac host galaxies. In this context, the U- 
and B-bands are especially important as they provide photometry 



shortward of the 4000 A break, without which only approximate 
average ages can be determined. Note also that dust reddening 
has little effect, since any reddening only would imply even bluer 
intrinsic colours. The SED modelling utilizes all the colour in- 
formation simultaneously by fitting the SED of evolution syn- 
thesis model to calibrated fluxes of the host galaxies at different 
wavelengths. Each wavelength band represents one data point in 
the SED of the object. Comparing fits made for different model 
spectra gives an estimate of the ages of the dominant stellar pop- 
ulation component s in the host g alaxy. Such SED modelling has 
been performed bv lJahnke et al.l (12004.) for low redshift RLQ and 
RQQ hosts . For the analysis we used the P EGASE2 evolution- 
ary model (iFioc & Rocca-Volmerangell997 ) and a single stellar 
population (SSP) model of single metallicity because the age- 
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Fig. 7. The U-B colour profiles for BL Lac hosts. 



metallicity degeneracy in models cannot be resolved with mul- 
ticolour data only. We used instantanous burst models which as- 
sume that the young stellar population is formed in a short burst 
with an IMF and evolves passively thereafter. 

Synthesized galaxy spectra composed from libraries of stel- 
lar spectra can be used to derive the colours of the galaxies over 
a large wavelength range from UV to NIR for different ages and 
metalhcities . For generating the synthetic spectra we used the 
LScald(ll990h initial mass function, solar metallicity and five dif- 
ferent model ages (0.1, 0.7, 2, 6 and 14 Gyr), to allow a direct 
comparison with the results of 'J ahnke et al.l (|2004|) . 

From the synthesized spectra we have derived the U-B, B-V, 
B-R, R-H and B-H colours for each model ages and compared 
those with the observed colours of each BL Lac using an iter- 
ative least squares fitting, with age as the only free parameter 
The results of the fitting are presented in Table |7] where the best 
fitting model is marked as 1 and the poorest fitting model as 3. 

The fitting results indicate that the best fitting models are 
generally consistent with a young/intermediate age stellar pop- 
ulation. Only for two BL Lacs (lES 0229-H200 and MRK 421), 
the best fitting model is older than 2 Gyr and is best consistent 
with the 6 Gyr model. Note that these objects also have the red- 
dest host galaxies in the sample {B-H ~ 4.4, see Table |5]l. The 
model fitting is thus consistent with the host galaxy colours, sug- 
gesting that the host galaxies of these objects are dominated by 
an old stellar population, with no evidence for a young popula- 
tion. Eight objects are best fit with the 2 Gyr model, while the 
remaining eight objects in the sample require an even younger 
model, 0.7 Gyr, and in the case of MS 1552.1-H2020, the 0.1 
Gyr model. Note that MS 1552.1+2020 deviates from the sam- 
ple also in the Kormendy relation (see above). 

These results are in excellent agreement with similarly 
young stellar populations found by Jahnke et al. (2004) for RLQ 
and RQQ host galaxies obtained from population synthesis 
model. They found the 2 Gyr model to best fit their data in 



Table 7. Colours of the single SSP fit compared to the colours 
of the BL Lac sample. 1 represents the best and 3 the worst fit. 



Name 


0.1 


0.7 


2.0 


6.0 


14.0 




Gyr 


Gyr 


Gyr 


Gyr 


Gyr 


lES 0229+200 






3 


1 


2 


PKS 0521-365 


3 


1 


2 






PKS 0548-322 


2 


1 


3 






IH 0706+591 




2 


1 


3 




MRK 421 






2 


1 


3 


MRK 180 


3 


1 


2 






lES 1212+078 






1 


2 


3 


MS 1229.2+6430 


2 


1 


3 






lES 1255+244 


2 


1 


3 






PG 1418+546 




2 


1 


3 




lES 1426+428 


3 


2 


1 






lES 1440+122 






1 


2 


3 


AP LIBRAE 


3 


1 


2 






MS 1552.1+2020 


1 


2 


3 






MRK 501 




2 


1 


3 




IZw 187 




2 


1 


3 




3c371 




2 


1 


3 




BL LACERTAE 


3 


1 


2 







12 out of their 19 objects, and an even younger population (0.7 
Gyr) in five objects, while an older population (6 Gyr) was pre- 
ferred in only two objects. Note that, as is the case for the quasar 
hosts (Jahnke et al. 2004), none of the elliptical BL Lac hosts 
are best modeled with the 14 Gyr model, i.e. a very old, evolved 
population as would be expected for early-type host galaxies. 
On the other hand, for both types of AGN hosts, there is also 
little evidence for massive ongoing starbursts with a signifi- 
cant very young population (age << 1 Gyr). Together these re- 
sults are in good agreement with the colour information (sec- 
tion 3.3) and strongly support the idea that the blue colours 
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Fig. 5. Upper: The B-R vs. 7? colour-magnitude diagram for BL 
Lac hos t galaxies (filled ci rcles), RLQ and RQQ hosts (open 
circles; Jahnke et al.l (12004 ). RGs (open squares; Govo ni et al.i 
(2000)) and inactive early-ty pe galaxies (asterisks; Pelet ier et al.l 
(il99a) ; iColbert et al.1 (l200lb ). Lower: The B-H vs. H colour- 
magnitude diagram for BL Lac host galaxie s and c omparison 
samples. Open circles are fr om iJahnke et alJ (12004 and aster- 
isks from iBoweretalJ (119921) . 



of the early-type host galaxies of moderately luminous AGN 
are caused by them having experienced a relatively recent SF 
episode. Further support for this conclusion is provided by spec- 
troscopic evidenc e for young/intermediate age populati ons in 
AGN hosts (e.g., Nol an et all 120011: iRaimann et all 1200 5). On 
the other hand, both th e predominantly late-type hosts of low 



luminosity AGN (e.g.. iKotilainen & WardI [19941: ISchade et all 



l2000t I Jahnke et alJl2004 and the elUptical hosts of very lumi 



nous AGN (e.g.. lDunlopetani2003h exhibit similar colours to 
their inactive counterparts. 

The blue colours and steep colour gradients of early-type 
AGN hosts found in an increasing number of studies are most 
likely caused by a young stellar population, and indicate a link 
between SF and the onset of the nuclear activity, both likely trig- 
gered by a tidal interaction or a minor or major merging event. 
However, the lack of obvious signs of interaction (companion 
galaxies, tidal tails, disturbed morphology) in the close environ- 
ment of the large majority of the host galaxies may require a 
significant time delay (at least hundreds of Myr) between the 
event with associated SF episodes and the start of the nuclear 
activity. Such a time delay is indeed predicted by simulations of 
gas reaching th e galaxy center after a merger/interaction (e.g., 
iLin et al.l[T988h . 

4. Conclusions 

We have presented Z?-band imaging of a sample of 18 low red- 
shift (z < 0.3) BL Lac objects for which their host galaxies 
were previously resolved in the optical R- and NIR //-bands. 
Subsamples consisting of nine and five BL Lacs were also im- 
aged in the U- and V-bands, respectively. All the BL Lac ob- 
jects were clearly resolved in the B- and V-bands while 8/9 
objects were resolved in the //-band. All the host galaxies are 
well described by an elliptical de Vaucouleurs model. These data 
are combined with previous optical and NIR data to study the 
optical-NIR colours and colour gradients of the host galaxies. 

BL Lac host galaxies are luminous (massive) and large ellip- 
tical galaxies having average B-band absolute magnitude Mb = 
-21.6 + 0.7 and average effective radius R^ - 7.6 + 3.2 kpc. 
They are indistinguishable from low redshift radio galaxies but 
brighter than inactive elliptical galaxies. The B-band Kormendy 
relation for BL Lac hosts is slightly steeper than that of RGs 
but does not deviate significantly from the Fundamental Plane 
of inactive elliptical galaxies. This indicates that BL Lac hosts 
are dynamically similar to normal ellipticals and that the active 
nuclear phase does not have any significant influence on the dy- 
namical structure of the galaxy. Thus it is possible that all mas- 
sive galaxies can experience an AGN phase. 

The distributions of the integrated blue/NIR colours (B- 
H = 3.5 + 0.5; B-R = 1.2 + 0.5) and colour gradients (A(B- 
R)/A(logr) - -0.14 + 0.75) of the BL Lac hosts are much wider 
than those for normal ellipticals with old stellar populations, and 
many BL Lac objects have bluer hosts and/or steeper colour 
gradients than those in normal ellipticals. The blue colours and 
steep colour gradients are most likely caused by a young stel- 
lar population, and indicate a link between SF caused by an 
interaction/merging event and the onset of the nuclear activity. 
Many targets have inverted (positive) colour gradients which in 
some cases, e.g. lES 1426+428, are most likely caused by ra- 
dial variations in the dust content of the host galaxy, but in some 
other cases may suggest star formation in the inner regions of 
the host galaxy. This result is corroborated by stellar population 
modelling, indicating young/intermediate age populations in the 
majority of the sample, in agreement with low redshift quasar 
hosts. However, the lack of obvious signs of interaction may re- 
quire a significant time delay between the event with associated 
SF episodes and the start of the nuclear activity. 

Future work in this area should address the correlation be- 
tween colour and statistics of companion galaxies and morpho- 
logical disturbances as an indicator for the interaction, for well 
defined samples of AGN hosts and inactive galaxies. In a forth- 
coming paper, we shall present NIR spectroscopy of BL Lac 
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hosts and RGs with blue colours, to analyze their stellar con- 
tent and SF properties, based on emission and absorption line 
diagnostics, in more detail than is affordable with imaging. 
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Fig. 8. The observed f/-band radial luminosity profiles (solid points with error bars) for each BL Lac, overlaid with the PSF model 
(dotted line), the de Vaucouleurs r'^"* model (long-dashed line) and the fitted PSF+host galaxy model profile (solid line). The X-axis 
is in arcsec and the Y-axis in mag arcsec 2. 
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Fig. 9. The observed fi-band radial luminosity profiles and model profiles. For explanation, see the caption of Fig. [8l 
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Fig. 10. The observed V-band radial luminosity profiles and model profiles. For explanation, see the caption of Fig. [8] 



